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I. INTRODUCTION 

Within the standard model (SM) of particle physics, 
CP violation in the quark sector is described by the 
Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing ma- 
trix. Physics beyond the SM may result in measured val- 
ues of observables that deviate from the values expected 
based on other CKM parameter measurements and the 
SM. 

The decay i?" — >■ Tr+Tr^Tr" [1] is well suited to the study 
of CP violation and has been previously explored by both 
the BaBAR [2] and Belle [3] Collaborations. Early studies 
of this mode were "quasi-two-body" (Q2B) analyses that 
treated each p resonance separately in the decays B^ -^ 



p"7r0(p" -^ TT+TT-) and S" 



p^TT 



=F(p± -^ 7r±7r"). How- 
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ever, as first noted by Snyder and Quinn [4j, a complete 
time-dependent Dalitz plot (DP) analysis is sensitive to 
the interference between the strong and weak amplitudes 
in the regions where the p"*", p~ , and p° resonances over- 
lap. This interference allows the unambiguous extraction 
of the strong and weak relative phases, and of the CP- 
violating parameter a = aT^gi-V^^V^l/ (V^jV*^)], where 
Vqqi are components of the CKM matrix. A precision 
measurement of a is of interest because it serves to fur- 
ther test the SM and constrain new physics that may 
contribute to loops in Feynman diagrams. 

In this paper, we present an update of an earlier BaBAR 
analysis. We use the full BaBAR dataset collected at the 
T(45') resonance, corresponding to an increase of 25% in 
the number of B meson decays, and include a number of 
improvements to both the reconstruction and selection 
procedures. Among these are improved charged-particle 
tracking and particle identification (PID), and a reop- 
timized multivariate discriminator, used both for event 
selection and as a variable in the final fit. 

Section [H] contains an introduction to the theory be- 
hind this analysis and the formalism used. We proceed 
to descriptions of the detector (Sec. HI), the datasets 
(Sec. IV), and the event selection procedures (Sec. [v|). 
This is followed by a presentation of the fitting proce- 
dure (Sec. VI I and of the systematic studies (Sec. VII I. 



Finally, we present the fit results (Sec. VIII I and a con- 



clusion (Sec. IX I. An overview of robustness studies is 
provided in an Appendix. 



II. THEORY OVERVIEW 
A. Time-Independent Probability Distribution 

The time-independent amplitudes for B'^ and B de- 
cays to 7r"'"7r~7r'' are given by 



As^ = f+T + fJA- + f„A\ 



(1) 



respectively, where A" and A with k e {+,—,0} are 
complex amplitudes associated with the /o"*", p~ , and 
p" resonances, respectively, and /„ = fnimjOi^) are de- 
fined in terms of modified relativistic Breit-Wigner reso- 
nances [S] modeling the three p resonances. The angle 0^ 
is the helicity angle for the resonance, defined as the an- 
gle between the 7r° (tt" ) momentum and the negative of 
the momentum of the recoiling 7r~ (7r+) for the /o+ {p~), 
and as the angle between the tt"*" momentum and the neg- 
ative of the momentum of the recoiling 7r° for the p" . All 
helicity angles are calculated in the p rest frame. In the 
fit, we include the p{770) as well as its radial excitation, 
the p(1450); therefore, each /^ is a sum of modified rela- 
tivistic Breit-Wigner resonances, F, for the p'^{770) and 
p'^(1450): 

fK{m,9^) ex i^p(770)(™,^K) + Ve"^'''^p(1450)(TO,6'«), 

(2) 

where Op' and (j)pi are the magnitude and phase of the 
p(1450) resonance relative to the p(770). We include sys- 
tematic uncertainties, described in Sec. |VII A} to account 
for possible contributions from the p(1700). 



B. Time-Dependent Probability Distribution 

Using the time-independent amplitudes A^j^ and A^^^, 
we can express the full time-dependent probability for a 

meson that is a B^ (-^3^) or B (A^^) at the time the 
other B meson decays, to decay to 7r"'"7r~7r'' as 

g-|At|/r^o 



i^r.(Ai)i^ 



4tro 



As.r + lAs.l 



T (l^a^P - \A3.n cos{AmdAt) 

'q 



±2 Im 



4o A* 



P 



sin(AmdAt) 



(3) 



where Tgo is the mean neutral B lifetime, Am^ is the 
mass difference between the heavy and light neutral 



B mass eigenstates, p and q are the complex param- 
eters in the definitions of the neutral mass eigenstates 
p\B^) ± q\B ), and At is the time difference between 
the decays of the fully reconstructed B meson {B^y^) and 
the B meson used to determine the B flavor (Stag)- In 
Eq. ([3]), as in the fit, we assume that the heavy and light 
mass eigenstates have the same lifetime, that there is no 
CP violation in B°B° mixing {\q/p\ = 1), and that CPT 
is conserved. 



C. Square Dalitz Plot Formalism 

While nonresonant phase-space decays uniformly pop- 
ulate the kinematically allowed region of a DP, signal 
pn events populate the boundaries of this region due to 
the low mass of the p resonances relative to the B mass. 
In particular, the interference regions of the signal DP, 
which provide sensitivity to the relative phases of the p 
resonances, are confined to small regions in the three cor- 
ners of the DP. In order to expand these regions of inter- 
est and avoid the use of bins of variable size, we perform 
a transformation of the DP that maps the kinematically 
allowed region onto a dimensionless unit square. The 
transformation is described by 

dm+dm^ -^ \deiJ\dm! dO' , (4) 

with the square Dalitz plot (SDP) coordinates 



1 

TT 
1 



Too 



- 1 



(5) 
(6) 



where m± is the invariant mass of the tt tt" system, toq is 
the invariant mass of the two charged pion candidates, ^o 
is the p° helicity angle defined earlier, m™^^ — mgo —m^o 
" — 2m^+ are the kinematic limits of the toq 



and 

mass, and J is the Jacobian of the transformation. 

determinant of the Jacobian is given by 



|detJ| =4|p;||p5|mo 



dmo d cos Oq 
dm! dO' 



where 



\/(^ 



ip;i 



pSI = V(^o*)'-<o, 



The 

(7) 

(8) 
(9) 



and the energies E'^ and -Eq of the tt"*" and 7r° are defined 
in the tt+tt" center-of-mass (CM) frame. Figure [l] shows 
an example of a standard DP (left) and its transformed 
SDP counterpart (right), plotted using simulated B^ — >• 
piT decays, where the three p resonances are assumed to 
have the same amplitude. 



1^ B% TiVn" (kin.) 
^ □ interference regs. J::-:?ias:; 
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FIG. 1: (color online) Nominal (left) and square (right) B° — >■ p7r Dalitz plots obtained from Monte Carlo generated events 
without detector simulation [2]- The amplitudes in Eq. ([I| are generated with values A^ = A~ = A" = 1 so that they interfere 
destructively for equal tttt masses. The hatched areas indicate the main overlap regions between the different p bands. The 
dashed lines in the square Dalitz plot correspond to m(p"'"'~''') — 1.5 GeV/c^. The middle plot depicts the Jacobian determinant 
of the transformation and shows the distribution in the square Dalitz plot for uniformly distributed events in the nominal Dalitz 
plot. 



D. U/I Formalism 

If one explicitly inserts Eq. (IT]) into Eq. (Is]), the full 

time-dependent amplitude for a B'^ or B meson to decay 
to 7r+7r^7r° can be written in terms of 



1^3x1' ±1-43. 



Ke[+-,a] 

+2 E (Re [/./;] f/±/'^ 

-Im [/./;] C/±,i-), 



(10) 



and 



Im 



-^3.^3^ 



K.e[+,-,o] 

+ E {^^if^f:]it 

fc<cre[+,-,0] 



(11) 



with 



u^ 


= lA" 


2± 


3'!', 




(12) 


rr±,Ro(Im) 
"-^Kcr 


= Re(Im) 


A'^A"* ± Ta" 


, (13) 


I. 


== Im 


A" A"* 


7 


(14) 




= Re 


A'" A"* - A" A''* 


; 


(15) 


rlin 


= Im 


a" A"* +A" A''* 




(16) 



The 27 real-valued U and / coefficients provide an al- 
ternative parameterization to tree and penguin ampli- 
tudes (as well as a) or to the amplitudes A" and A |6|. 



The U and / parameters can also be directly related 
to the Q2B C and S parameters often used in CP- 
violation analyses [7], where C parameterizes direct CP 
violation, and 5 parameterizes mixing-induced CP vio- 
lation (involving the angle a in this analysis). The re- 
lated parameter AC describes the asymmetry between 

the rates r(B" -^ p+n-) + r(B" -^ p-'K+) and r(B" -^ 
p^TT^) +V{B — > p+vr^), while A5 is related to the 
strong phase difference between the different amplitudes 
describing the decay B^ — >■ p-K. The U and / parame- 
ters are related to the C and S parameters through the 
relations 



c-^^.c-^^- 



ui 



ui 






2/_ 



and 



where 



A -^+-^- 



uX + ur 



(17) 



(18) 



C = {C+ +C-)/2, 


(19) 


AC = {C+ -C-)/2, 


(20) 


S={S++S-)/2, 


(21) 


AS= {S+ -S-)/2. 


(22) 



Note that while C, S, AC, and AS do not depend on 
interference effects between the p resonances, the U and 
/ parameter formalism accounts for these features and 
is thus appropriate for a full DP analysis. While some 
degree of physical intuition is lost when using the U and 
/ parameters instead of the standard complex amplitudes 



and phases, there are several practical motivating factors 
for their adoption in the fit: 

• Whereas there is a two-fold ambiguity in a (a ver- 
sus 90° — a), there is a unique solution in a fit to 
the U and / parameters, which encompasses both 
solutions for a). 

• The U and / parameterization results in uncertain- 
ties that are more Gaussian than those in a stan- 
dard amplitude and phase fit. 

• It is simpler to average U and / fit results from 
different measurements or experiments that publish 
the full covariance matrix. 

For physical solutions, there are constraints between 
the U and / parameters. In the case of three p reso- 
nances (p+, p^ , and p°), a fit to the complex tree and 
penguin amplitudes as well as the weak phase a involves 
11 unknown parameters, which reduce to 10 parameters 
when the arbitrary global phase is removed. A U and 
/ fit is equivalent to such a fit, but involves many more 
parameters. However, when the p'^n'^ amplitude is small, 
as is observed in nature, the values of 11 of the 27 U 
and / parameters become unimportant. Due to the high 
degree of correlation between the various U and / fit pa- 
rameters, the 27 parameters actually represent only 12 
independent parameters. Neglecting the arbitrary phase 
and the overall normalization, this reduces to 10, and 
once isospin relations are taken into account, only 9 in- 
dependent parameters remain. 

Because the U and / formalism is used in the final fit 
without any constraints on the parameters (aside from 
fixing U^ = 1 to set the overall normalization), it is 
possible for the free parameters to take on unphysical 
values that do not correspond to any physical set of p 
amplitudes. The final fit values from the 2007 BaBar 
analysis |5j are one such unphysical set. We determined 
that no biases are introduced due to the fitted values of 
the parameters being unphysical. 



III. THE BABAR DETECTOR AND 
EXPERIMENT 

The data used in this analysis were collected with the 
BaBar detector at the PEP-II asymmetric-energy e^e" 
storage ring at SLAG. A detailed description of the BaBar. 
detector is presented in Ref. [S]. The tracking system 
used for track and vertex reconstruction has two compo- 
nents: a 5-layer silicon vertex tracker and a drift cham- 
ber, both operating within a 1.5 T magnetic field gen- 
erated by a superconducting solenoidal magnet. A de- 
tector of internally reflected Gherenkov light associates 
Gherenkov photons with tracks for particle identifica- 
tion. The energies of photons and electrons are deter- 
mined from the measured light produced in electromag- 
netic showers inside a GsI(Tl) crystal electromagnetic 



calorimeter. Muon candidates are identified with the 
use of the instrumented flux return of the solenoid. The 
flux return instrumentation initially consisted of resistive 
plate chambers and was later modified to consist of a 
mixture of resistive plate chambers and limited streamer 
tubes. 



IV. DATA SAMPLE AND MC SIMULATION 



A. Data Samples 



For the final fit, we use the full "on-resonance" BaBar 
dataset of 431.0 fb^^ collected at the T(4S') resonance 
energy {y/s = 10.58 GeV/c^). When optimizing back- 
ground suppression criteria, 44.6 fb~ of "off-resonance" 
data, collected 40MeV below the T(4S') resonance, are 
used to model "continuum" e+e^ ^ QQ {q = u,d,s,c) 
background. 



B. Monte Carlo Samples 

Event simulation based on the Monte Garlo (MG) 
method is used to evaluate backgrounds and to determine 
signal-event reconstruction efficiencies. For all MG sam- 
ples, detector response is accounted for using the GEANT4 
package [9 in a full BaBar, detector simulation. 

B° — 7> p+TT^ and B'^ — > p^Tr^ signal decays are simu- 
lated in separate MG datasets, but the B" -^ p^Tr^ sam- 
ple is not used to determine background selection criteria 
or to model signal distributions since the nominal branch- 
ing fraction for B° — >■ p^n^ is 11.5±3.1 times larger than 
the branching fraction for B° —>■ p^ir'^ [lO'. 

B-decay backgrounds are modeled using MG samples 
consisting of B decays to specific final states as well as 
"generic" MG samples consisting of charged and neutral 
B decays to unconstrained final states. In the generic 
MG, dominant branching fractions are fixed to the results 
of recent measurements [H] . Due to the uncertainty on 
branching fractions for charmless decays, all charmless 
events are removed from the generic MG and charmless 
modes of interest are explicitly included among the B 
background samples consisting of decays to specific final 
states. We use a total of 24 B-decay MG samples corre- 
sponding to 29 different final states (Table IT]) as well as 
MG samples of generic charged and neutral B decays. 

The expected number of events for each charmless B 
background is calculated according to 



— "^nsB^^nwdc^, 



(23) 



where ubb is the number of produced BB pairs, B is the 
branching fraction [TT] (approximately 1/2) for an T(4S') 
to decay to a charged or neutral BB pair (whichever 
is appropriate for the mode in question), yBmodc is the 
branching fraction for the B decay mode, and e is the effi- 
ciency for reconstructing events in the mode, determined 



TABLE I: B background decay modes included in the final fit. Modes generated taking into account interference effects in the 
Dalitz plot are indicated by the "Dalitz" label. Longitudinal polarization is indicated by "[longitudinal]". 



Class 


Decay mode 




B [10-"] 


# of expected events 





B+ 


— > p^p" [longitudinal] 




24.0 ±1.9 


129 ± 10 





B+ 


^a+(^(p7r) + y 




26 ±7 


53 ±14 





B+ 


-> a?(-> p±7rT)7r+ 




20 ±6 


37 ±11 


1 


B+ 


-^ 7r+7^^'(^ TT+TT-) 




7.99 ±0.35 


6.96 ±0.31 


2 


B+ 


-> A'+tt+tt" Dalitz 




51.0 ±2.9 


34.8 ±2.0 


2 


B+ 


— 5> TT+TT+TT^ Dalitz 




16.2 ±1.5 


203 ± 19 


3 


B+ 


->^«p+ 




10.9 ±1.4 


120 ±15 


3 


B+ 


^ ^+7^^'(^ ,r'V") 




3.54 ±0.15 


24.2 ±1.1 


4 


B+ 


^ ^+^" 




5.7 ±0.5 


38.6 ±3.4 


4 


B+ 


^ T^+tt" 




12.9 ±0.6 


18.6 ±0.9 


5 


B" 


-> K^+Tv' Dalitz 




44.8 ±2.6 


15.7 ±0.9 


6 


B" 


— > p^p^ [longitudinal] 




24.2 ±3.1 


122 ±16 


6 


B° 


— > aj^7r+ 




33 ±5 


61 ±9 


6 


B° 


^ aV" 




11.0 ±1.7 


22.8 ±3.5 


7 


B" 


^K+TT- 




19.4 ±0.6 


21.6±0.7 


8 


5° 


-^ K+TT-n" Dalitz 




35.9 ±2.6 


398 ± 29 


9 


5° 


^ 7/s:*"(892)(^ K+Tv- 


) 


40.1 ±2.0 


31.8±1.6 


9 


B" 


-^-fK'-^\U30){^ K+TT 


") 


12.4 ±2.4 


3.2 ±0.6 


9 


B" 


^V(^p"7y' 




0.35 ±0.18 


4.2 ±2.1 


10 


B" 


^7r"A:^^(^7r+7r-) 




3.39 ±0.21 


21.8 ±1.4 


11 


B" 


^B-(^7r-7r")7r+ 




3.35 ±0.27 


399 ± 32 


12 


B" 


-^d"{^K+7t-, K+7V- 


-vv« 


46.7 ±4.5 


124 ±12 


13 


B" 


-^d"{^tv+tt-)tv'' 




0.367 ±0.034 


48 ±5 


14 


B" 


^J/V(->e+e-,^V") 


^« 


2.09 ±0.19 


153 ±14 


15 


B° 


— > neutral generic 6 — >• c 


decays 




466 ± 14 


16 


B+ 


— >• charged generic b — > 


c decays 




921 ±21 




Total 






3478 ± 65 



from MC. The factor of 2 is included because either of 
the B mesons in a given event may decay to the mode of 
interest. 

In the case of the charged and neutral generic B back- 
grounds, the number of events expected for each mode 



nuc- 



ubbB 



(24) 



(/gen 



V. EVENT SELECTION AND BACKGROUND 
SUPPRESSION 

A. Event Preselection 

The kinematics of B meson decays that are fully recon- 
structed at BaBar can be characterized by two variables: 
Toes and AE. The beam-energy-substituted mass to.es 
is the invariant mass of the reconstructed B candidate 
calculated under the assumption that its energy in the 
e^e" CM frame is half the total beam energy. We define 



where rigon is the number of generated charged (neutral) 
MC events and umc is the number of charged (neutral) 
generic MC events that remain after all selection criteria 
have been applied. 

An additional simulated dataset, which is used for val- 
idation studies, consists of DP-parameterized B^ -^ pn 
decays. This dataset is used to verify fiavor-tagging con- 
ventions. 



TOES 



[{s/2+p,-pB)/E,f-\pB\' 



(25) 



where 



/s is the total beam energy in the e^ 



fi- CM frame, 
(Ei^Pi) is the four-momentum of the e^ e' system in the 
laboratory frame, and pb is the _B-candidate momentum 
in the laboratory frame. The second kinematic variable 
is defined by 



AE^ E 



1 



B 



(26) 



10 



where E'^ is the measured energy of the B candidate in 
the e+e~ CM frame. 

Pairs of oppositely charged tracks are combined with 
7r° — > 77 candidates to construct B^ — > 7r+7r~7r° candi- 
dates. During a preselection stage, we require Eg to lie 
between 4.99 and 5.59 GeV. For the charged tracks cor- 
responding to the TT candidates, we require a maximum 
momentum of 10 GeV/c, a minimum of 12 hits in the 
drift chamber, a maximum distance-of-closest-approach 
(DOCA) relative to the beamspot center in the x-y plane 
of 1.5 cm, and a DOCA along the z axis between —10 
and -1-10 cm. We require the two photon candidates that 
form the 7r° candidate to have an energy in the laboratory 
frame between 30 MeV and 10 GeV, and lateral moments 
in the electromagnetic calorimeter less than 0.8. We re- 
quire a 7r° — T' 77 mass between 100 and 160 MeV/ c^, and 
a total 77 laboratory energy greater than 200 MeV. 



B. Primary Selection Criteria 

Following preselection, we require toes to lie between 
5.200 and 5.288 GeV/ c^. For the charged tracks corre- 
sponding to the TT^ candidates, we require a minimum 
transverse momentum of 0.1 GeV/c. The lateral moments 
of the two photons from the 7r° candidate are each re- 
quired to lie between 0.01 and 0.60, while the laboratory- 
frame energies in the electromagnetic calorimeter are re- 
quired to exceed 50 MeV. The mass of the 7r° candidate 
(TO-y-y) is required to satisfy 110 < m^o < 160 MeV/ c^. 

We calculate At by measuring the distance along the 
beam axis between the B:^.^ and Stag decay vertices and 
using the boost {/3j « 0.56) of the e+e~ system. The 
time difference At and its estimated uncertainty are re- 
quired to satisfy \At\ < 20 ps and a{At) < 2.5 ps. 

The average number of B candidates (measured as the 
total number of candidates that pass all the preceding 
selection criteria divided by the total number of events 
with at least one candidate) is 1.45. To retain only one 
-B-decay candidate in each event, we select the candidate 
that has a tt" — )■ 77 candidate mass closest to the world 
average value of the 7r° mass. In the case of multiple B 
candidates reconstructed with the same tt" candidate, we 
select one B candidate arbitrarily. A tighter requirement 
of 5.272 < Toes < 5.288 GeV/ c^ is applied after selecting 
a single candidate in each event. 

We use a lower toes sideband of on-resonance data, 
5.215 < Toes < 5.270 GeV/ c^, as well as off-resonance 
data, to model the distribution of continuum events. 



AE±{m, 



c± ~ (c± T c)(to^+^- /m^^l. )2, to^4 



is the measured tt+tt^ mass, and the parameter to™+^_ 
is fixed at 5.0 GeV/c^ (corresponding roughly to the max- 
imum observed value of to^+^-). The c and c± param- 
eters are calculated from fits to B^ — >■ p'^iT~ signal MC. 
The dimensionless quantity AE' serves to reduce the de- 
gree of correlation with m^+j^- and is included as an 
input variable in the final fit. 

To calculate the three c parameters, signal MC events 
are divided into seven equal-sized bins in m^+T^- from 
to 5.143 GeV/ c^. In each mass bin, the peak in AE is fit 
with the sum of two Gaussians, and the mean (/i) and 
width (a) of the narrower Gaussian are extracted. From 
these parameters, two sets of datapoints are constructed: 
one consisting of the values oi jj. + 3a for each mass bin, 
and the other consisting of the values of /i — 3(T for each 
mass bin. The first set is fit with the quadratic function 
corresponding to AE^{mT^+^~), which yields values for 
c+ and c. A similar fit is performed on the second set 
using the function A£'_(TO7r+7r-)j which yields values for 
c_ and c. The c value used in the final transformation is 
obtained by averaging the c values from the separate fits. 
For the final calculation of AE' , we use the parameter 
values c = 0.0792 GeV, c_ = -0.1433 GeV, and c+ = 
0.1093 GeV. Candidates are selected if AE' is between 
— 1 and -fl, which is roughly a ±3cr criterion on AE. 



D. 



Particle Identification Selection Criteria 



In the previous BaBar. analysis [2], the charged pion 
candidates were required to be inconsistent with muon, 
kaon, proton, and electron hypotheses. Improvements 
in the BaBar, reconstruction software now provide de- 
creased false-positive rates for a given signal efRciency. 
We apply a PID selection criterion corresponding to a 
signal-to-square-root-of-background (S/^/B) ratio (cal- 
culated using off-resonance data to model background, 
and correctly reconstructed B° — >■ p^n^ signal MC to 
model signal) of 1.26 (scaled relative to all previous se- 
lection criteria). For comparison, application of the PID 
criteria of the previous analysis achieves a scaled S/^/B 
ratio of 1.14. The new PID criterion selects 54.6% of 
off-resonance data and 93.2% of correctly reconstructed 
signal MC events relative to the previous selection crite- 
ria. 



E. 



Multivariate Discriminator Selection 



C. Transformed AE' Definition and Selection 
Criterion 

Because the width of the AE distribution is highly 
correlated with the tt+tt" mass and hence varies across 
the DP, we introduce the dimensionless transformed vari- 
able AE' = (2A£; - AE+ - AE^)/{AE+ - AE^ ) , where 



Continuum qq events are the primary source of back- 
ground in this analysis. To improve discrimination be- 
tween continuum background and signal events, we train 
a neural network (NN) on off-resonance data and signal 
MC using the "MLP" (multi-layer perceptron) NN im- 
plementation in the TMVA software package pj]. The 
NN is trained using 10,000 events representative of sig- 
nal, 7,500 events representative of background, and 200 
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training iterations. Validation is performed using an in- 
dependent sample of 20,000 signal and 7,500 background 
events. The NN is configured to use two hidden layers 
with 6 and 5 nodes, respectively. The NN uses four event- 
shape variables that help distinguish between the roughly 
isotropic shape of B decays and the more jet-like shape 
of continuum events. The training variables include the 
Legendre moment LO of the rest of the event, defined as 
the sum of the magnitudes of the momenta of all charged 
particle candidates not belonging to the reconstructed B 
candidate; the ratio of the Legendre moment L2 to the 
Legendre moment LO, where L2 is defined as the sum 
over all tracks and neutral clusters not belonging to the 
B candidate oi p ■ (3 cos^ 6 — l)/2 where p is the magni- 
tude of the momentum of each track or cluster and is 
the angle between the B thrust axis and the momentum 
corresponding to the track or cluster; the cosine of the 
angle between the B candidate momentum and the beam 
axis; and the cosine of the angle between the thrust axis 
of the B candidate and the beam axis. All of these input 
variables are calculated in the CM frame. 

While optimizing the NN discriminator, we studied 
whether some improvement in performance might be 
achieved by training the discriminator separately in each 
of seven B-flavor tagging categories (each of which has a 
different flavor tagging purity), or applying NN selection 
criteria separately in each of the _B-flavor tagging cate- 
gories. These approaches were found to yield negligible 
improvements in signal efficiency for a given degree of 
background rejection. 

The final performance of the NN discriminator is 
shown in Fig. [2] while Fig. |3] shows the separation 
achieved between signal and background samples dur- 
ing training. The NN is used to apply a loose selection 
criterion that retains 75% of the signal events remaining 
after all previous selection criteria. The NN output is 
also used as an input variable in the final fit. 



Selection Performance 
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FIG. 2: (color online) S'^ — >■ p tt^ continuum background 
rejection versus signal efficiency for final neural network im- 
plementation, with continuum background represented by off- 
resonance data, and signal represented by signal MC. Points 
on the curve correspond to different requirements on the neu- 
ral network output. Our selection criterion retains 75% of 
signal events. 
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6 - 



Background 




Table ITT] provides a summary of the signal efficiencies 
for each step of the selection process for simulated B'^ — > 
p+7r~, B^ — >• p°7r°, and nonresonant (three-body phase 
space) B'^ — > 7r+7r^7r° events. 



FIG. 3: (color online) NN output distributions for signal 
B" — >■ p^TT^ (solid blue) and continuum background (hatched 
red) with an arbitrary vertical scale. Continuum background 
is represented by off-resonance data while signal is represented 
by signal MC. Our selection criterion requires NN > 0.58. 



VI. MAXIMUM LIKELIHOOD FIT 



We perform an unbinned extended maximum likeli- 
hood fit to the selected events in order to extract the 
event yields and the U and / parameters. The input vari- 
ables are tties, AE' , the NN output, and the three timc- 
dependent-SDP variables m', 9' , and At. We also use 
(TAt (the per-event uncertainty on At) as a scale factor 
in the signal At resolution function. The likelihood func- 
tion used in the fit consists of separate components for 
signal, continuum background, charged B backgrounds, 



and neutral B backgrounds. The signal PDF is subdi- 
vided into components describing correctly reconstructed 
or "truth-matched" (TM), and misreconstructed or self- 
cross-feed (SCF) candidates. Truth-matched events are 
identified using MC generator information and contain 
correctly reconstructed 7r+, tt", and 7r° candidates, two 
of which were produced in a p decay. It is also required 
that the p and the remaining n come from the same B- 
meson parent. For nonresonant MC, we instead require 
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TABLE II: Selection efficiencies, ex , relative to previous cri- 
teria for simulated B" — >■ X events {X = p'^iT~, p^tt", non- 
resonant 7r+7r~7r°). Statistical uncertainties on all efficiencies 
are less than 0.001. 



Cuts 


MC 






Preselection 


0.650 


0.582 


0.547 


PT >0.1GeV 


0.994 


0.986 


0.999 


jAii < 20 ps & a{At) < 2.5 ps 


0.966 


0.957 


0.960 


Photon energy and lateral mom. 


0.928 


0.960 


0.943 


110 < m(77) < 160MeV/c^ 


0.982 


0.983 


0.983 


PID 


0.920 


0.936 


0.928 


5.200 < jTiES < 5.288 GeV/c^ 


0.994 


0.995 


0.996 


-1 < AE' < 1 


0.819 


0.863 


0.871 


NN > 0.58 


0.749 


0.768 


0.762 


5.272 < TTiES < 5.288 GeV/c^ 


0.798 


0.826 


0.937 


Cumulative efficiency 


0.255 


0.265 


0.279 



background components, namely 



where 



/sCF^37r-SCF,i 






B+Ms-jl' B 



3 = 1 



)n 



(28) 



• N'i^ is the total number of tt'^-k 7r° signal events 
in the data sample (both TM and SCF); 



that all three pions come directly from the same B par- 
ent with no intermediate p resonance. All signal events 



that are not TM are SCF. In B" ^ 



P^t:^ 



MC, 22% of 



reconstructed B candidates are SCF while in B^ 
MC, 14% are SCF. 



pOvrO 



The signal and _B-background components are further 
divided by B-fLavor tagging category and the flavor of the 
tag-side B. The Stag flavor is estimated using multivari- 
ate discriminator techniques and is classified as belonging 
to one of seven flavor-tagging categories corresponding to 
different degrees of probability that the flavor has been 
correctly determined [13j . 

Additionally, prior to the flnal fit, a transformation 
is applied to the NN variable. This has the effect of 
broadening the peak in the signal NN distribution, and 
transforming the continuum NN distribution to make fit- 
ting with the smoothing algorithm (see Sec. VI B) more 
effective. The transformation is defined by 



NN„ew - 1 



arccos(NNoid + /3), 



(27) 



where /3 = 0.029 is an offset approximately equal to 1 
minus the maximum of the original NN distribution. 



• /|^ is the fraction of signal events (TM and SCF) 
in fiavor tagging category c; 

• /sCF i^ the fraction of misreconstructed signal 
events (SCF) in tagging category c; 

• ^37T-TM,i ^'^^ ^3ir-scF.j ^'^^ ^^^ products of PDFs 
of the discriminating variables used in tagging cat- 
egory c for TM and SCF events, respectively; 

• N^-^ is the number of continuum events in flavor 
tagging category c; 

• ^tag.i is the tag flavor of the event where we use 
the convention that qta.g,i — 1 for -Btag = B'^ and 



9tag,; 



-1 for a B, 



tag 



#; 



A 



qq,taf. 



is the flavor tag asymmetry, parameterizing 
possible charge asymmetry in continuum events; 

• Vgq j is the continuum PDF for tagging category c; 

• ^cfass (-^cfass) ^^ ^^^ number of charged (neutral) 
i?-background classes included in the fit (including 
generic modes); see Table llj 

• Ng+j {Nqoj) is the number of expected events in 
the charged (neutral) B-background class j; 

• fs+j Wb°j) i^ ^^^ fraction of charged (neutral) B- 
background events of class j that are in flavor tag- 
ging category c; 



• ^s+,tag,j is the flavor tag asymmetry of charged 
i?-background class j; 



Likelihood Function 



• Vg+ j ■ is the i?+-background PDF for tagging cat- 
egory c and class j\ 



The PDF Vf for an event i in tagging category c is 
the sum of the probability densities of all the signal and 



• Vgo j is the B^-background PDF for tagging cat- 
egory c and class j. 
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The PDFs Vx are the product of PDFs for the dis- 
criminating variables xi = toes, X2 = AE', x^ — 
NNjiewj and the three time-dependent-SDP variables 
X4 = {m\e',At}: 



T'hu)^I['^hu)i^k). 



(29) 



A,-l 



The extended likelihood function including all tagging 
categories is given by 



N" 






(30) 



where N is the total number of events expected in tag- 
ging category c and iV^ is the observed number. 



B. Signal Parameterization 

The effect of experimental resolution for Ai in signal 
events is taken into account by convolving the PDF de- 
scribing the true Ai distribution with a sum of three 
Gaussians. 

The triple-Gaussian resolution function is constructed 
using a narrow "core" Gaussian, a slightly wider "tail" 
Gaussian, and a very wide "outlier" Gaussian. In the 
final fit to on-resonance data, all signal Ai parameters are 
fixed to values obtained from fits to fully reconstructed 
B decays. In the fits to fully reconstructed B decays, 
the mean and width of the outlier Gaussian are fixed 
to and 8 ps, respectively. Similarly, the width of the 
tail Gaussian is fixed to 3 ps, but its mean is allowed 
to vary. Finally, both the mean and width of the core 
Gaussian are allowed to vary in the fit, and the means of 
the core Gaussian are allowed to take on different values 
in each i?-fiavor tagging category. The means and widths 
of the core and tail Gaussians are scaled by the per-event 
uncertainty on Ai (cta*): 



'7^sig(At, CTAt) = 



(1 ^ /tail ^ /out) 

xG{At; S(.o„crAt, s^Tor^o-At) 

+ /tailG(Ai, Sja^iiCTAt; S^j^iiCTAt) 

-h/outG(Ai;Som,crout), (31) 



where G{x; xq, cr) is a Gaussian with mean xq and width 
a. See Fig. H (upper plot) for the distribution of At in 
signal MC events. 

The AE' signal distribution is modeled using the sum 
of two Gaussians where all five free parameters depend 
linearly on rn^(7r"'"7r~) in order to account for residual 
dependence on the mass. The to.es signal distribution is 
parameterized using a bifurcated Crystal Ball function, 
which is composed of a one-sided Gaussian and a Crystal 



Ball function: 

(j^-(x-mf/2sl fp^ (2; - to) > 0, 

fix) = { Ce-(^-™)'/2''i for -A< ^ < 0, 

Ci\)''e-^ {-^ - A - ^)-'' for ^ < -A. 

(32) 



The parameters describing both these lineshapes are ex- 
tracted from fits to signal MC events. In the final fit, the 
parameter m is free to vary for to,es while the core Gaus- 
sian mean and width, and the slope of the mean (i.e., the 
dependence of the mean on TO,^(7r+7r^)) are free to vary 
for AE'. 

The NNncw signal distribution is modeled by nonpara- 
metric histograms generated by smoothing the NNnew 
distribution in signal MC. 

The AE' distribution for SCF signal candidates is de- 
scribed using a single Gaussian with mean and width 
fixed to the values extracted from a fit to SCF signal 
MC. 

The distributions of toes and NNncw for SCF signal 
events are modeled with nonparametric histograms gen- 
erated by smoothing the appropriate one-dimensional 
distributions in SCF signal MC. The nonparametric 
PDFs for the NN„ew distributions in both TM and SCF 
signal are generated separately for each flavor-tagging 
category. 

While the SDP distribution for TM signal events is pa- 
rameterized by the full time-dependent decay probability, 
the SCF distribution is parameterized by modifying this 
distribution. Using signal MC, we create a binned map 
in the SDP that contains the probability for an event 
generated in each bin of the SDP to be reconstructed in 
the same bin, or each of the other bins. This map is con- 
volved with the time-dependent decay PDF to generate 
the SCF signal PDF in the SDP. 



C. Continuum Background Parameterization 

The At distribution for continuum background is mod- 
eled using a sum of three Gaussians. The parameters 
of this triple Gaussian are obtained from a fit to off- 
resonance data in which all parameters are allowed to 
vary. See Fig. [4] (lower plot) for the distribution of At in 
off-resonance data. 

The AE' distribution for continuum background is 
modeled using a second-order polynomial with param- 
eter values extracted from a fit to the to,es sideband in 
on-resonance data. 

The distribution of toes for continuum background is 
modeled using an Argus [14J function with shape and 
endpoint parameters that arc allowed to vary in the final 
fit. 

To account for residual dependence on DP position, 
the distribution of NNnow for continuum background is 
modeled using a second-order polynomial function in 
which each coefficient depends linearly on TO-dist, the 
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D. B-Background Parameterization 

The functional form of the At resolution functions for 
the B backgrounds is the same as that for signal. Pa- 
rameter values are obtained from separate fits to fully 
simulated MC data samples representative of each B- 
background class. 

Simulated samples for each of the ^-background 
classes are used to generate nonparametric PDFs for use 
in the final fit. One-dimensional PDFs are used for AE' , 
Toes , and NNnow , without any splitting by flavor-tagging 
category. Two-dimensional SDP PDFs are generated for 
each i3-background class and each B-flavor tag within 
that class. 



E. Dalitz-Plot-Dependent Selection Efficiency 
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Selection efficiencies across the SDP are calculated 
from a combination of all available nonresonant (_B" ^■ 
7r+7r-7r°), B° -^ p±7r=F, and B° -^ p°n° MC samples. We 
divide the SDP into a 40 by 40 grid and, for each bin, 
calculate both the fraction of events generated in that bin 
that are correctly reconstructed (TM), and the fraction 
of events generated in that bin that are misreconstructed 
(SCF). From this, we generate tables of efficiencies that 
are used as inputs to the fit. Histograms of the TM and 
SCF selection efficiencies are provided in Fig. [5] 



VII. SYSTEMATIC STUDIES 



A. Tlie Effect of the p(1700) 



FIG. 4: Distribution of At in simulated B" -> p tv^ signal 
events (top) and in off-resonance data (bottom). 



minimum invariant mass of any tttt combination in the 
S" ^- 7r+7r^7r° candidate, which acts as a measure of the 
distance from the edge of the DP. In addition, the poly- 
nomial is multiplied by (1 — NNnow )° where a is a linear 
function of TOdist • AH the polynomial parameters are free 
to vary in the final fit. 

The two-dimensional SDP distribution for continuum 
events is obtained by applying Gaussian kernel smooth- 
ing algorithms to the SDP distribution for on-resonance 
Toes sideband data and generating a two-dimensional his- 
togram from the resulting PDF, which serves as a non- 
parametric PDF in the fit. Histogram PDFs are gen- 
erated separately for each fiavor-tagging category and, 
within each category, for both S-fiavor tags. Bins in 
these histograms are mirrored across the 9' = 0.5 axis so 
that the distributions are symmetric in ^'. A number of 
parameters are allowed to vary in the fit to allow for an 
asymmetry. 



We include the p(1450) in the final fit with an assump- 
tion that the relative magnitudes and phases between the 
three p(1450) resonances are the same as for the p(770). 
Whereas there is reasonable motivation for this assump- 
tion in the case of the /9(1450) since the p{770) and 
/9(1450) have the same quantum numbers, the p(1700) 
does not share these quantum numbers (£ = 2 instead of 
0). Since the p(1700) is not expected to provide a large 
contribution, we exclude the p(1700) from the fit and 
associate a systematic uncertainty with this omission. 

Naively, one might calculate this systematic uncer- 
tainty by fitting the full dataset with and without the 
/9(1700) and calculating a single change in the U and / 
parameters, but this approach suffers from statistical un- 
certainties due to the fact that there is only one dataset 
available. Nonetheless, it is still useful to calculate a co- 
variance matrix using this approach in order to estimate 
the magnitude of the systematic uncertainty. As a first 
order assessment of the systematic uncertainty, we calcu- 
late a covariance matrix with elements given by 



a. 



d.6i 



(33) 



where Si is the difference between the value of the zth 
U or I parameter in fits with and without the p(1700). 
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TABLE III: Magnitude of changes in U and / parameters be- 
tween fits with and without the p(1700), the ratio of these 
magnitudes relative to statistical uncertainties from the fit 
without the p(1700), the bootstrapped estimate of the uncer- 
tainty on changes in U and I parameters between fits with 
and without the p(1700), and the ratio of the mean change in 
the U and / parameters across all bootstrapped fits to their 
bootstrap-estimated uncertainties. 



FIG. 5: (color online) Square Dalitz plot map of TM (top) 
and SCF (bottom) selection efficiencies. The color scale on 
the TM plot has a maximum of 0.60 while that on the SCF 
plot has a maximum of 0.075. 



Parameter 

To 

/_ 

rim 
-'-0 
rRc 

-f-o 

1+ 

rim 
-'+0 

rRc 
-'+0 

rini 

-'+- 

rRc 

-' + - 

Uo 

rr-,Rc 
TT+,Re 

uz 

u+ 

'^ + 

TT-,RC 
'^ + 

rr-t-.Im 
rr-^.Rc 

^-,Rc 
(y+Jm 

U+,R. 



\AU\ 

0.017 

0.015 

0.40 

0.39 

0.0013 

0.024 

0.18 

1.34 

0.68 

0.012 

0.006 

0.31 

0.31 

0.14 

0.25 

0.00041 

0.015 

0.18 

0.20 

0.039 

0.005 

0.6 

0.5 

0.038 

0.21 

0.021 



AU\/ast.t 




(At/) /aiV 


0.46 


0.020 


-0.36 


0.24 


0.022 


0.41 


0.92 


0.42 


-0.70 


0.67 


0.7 


0.28 


0.02 


0.023 


0.18 


0.07 


0.42 


-0.11 


0.40 


0.7 


0.26 


1.81 


0.9 


-0.95 


0.90 


0.9 


-0.59 


0.23 


0.029 


0.29 


0.21 


0.017 


0.54 


0.70 


0.5 


0.76 


0.89 


0.34 


0.49 


0.70 


0.21 


-0.76 


1.47 


0.17 


-0.77 


0.00 


0.034 


0.19 


0.21 


0.05 


0.15 


0.61 


0.44 


-0.10 


0.62 


0.36 


0.50 


0.24 


0.16 


0.51 


0.03 


0.17 


0.05 


1.23 


0.8 


-0.27 


1.10 


0.7 


0.15 


0.15 


0.25 


-0.22 


0.84 


0.28 


-0.21 


0.23 


0.06 


-0.46 



When the three /9(1700) resonances are included in the 
fit, their magnitudes and phases are all allowed to vary 
independently. The square roots of the diagonal elements 
of the matrix (in other words, the absolute value of the 
change in each U or I parameter) are given in Table III 
along with the ratio relative to the statistical uncertainty 
on each parameter in the fit without the p(1700). The 
ratios are generally less than 1; only 4 of the 26 ratios are 
greater than 1, indicating that the changes in the U and 
I parameters resulting from the inclusion of a p(1700) in 
the Dalitz model are mostly small. 

To assess the uncertainties on the changes in the U 
and / parameters, we employ the bootstrap technique 
introduced by Efron [18]. This approach allows us to 
calculate a covariance matrix associated with the uncer- 
tainty on the U and / parameters extracted from the fits 
to the full dataset by fitting a large number of datasets 
generated by sampling events with replacement from the 
full dataset. Each of the approximately 1000 resampled 
datasets generated in this manner has the same num- 



ber of events as the original dataset and is fit with and 
without the p(1700) component. For each pair of fits, 
the change in the U and / parameters is calculated and 
a covariance matrix is generated by determining the co- 
variances of the changes in the U and / parameters in 
all these pairs of fits. As demonstrated by Efron, the co- 
variance matrix obtained in this manner is an estimate 
of the covariance matrix associated with the uncertainty 
on the changes in the U and / parameters. Therefore, 
the square roots of the diagonal elements of this matrix 
are estimates of the uncertainty on the change in each 
U and / parameter. These uncertainties, as well as the 
ratio of the mean change of each U and / across all the 
bootstrapped fits to its estimated uncertainty (from the 
bootstrap), are given in the last column in Table III As 
the ratios in this table demonstrate, the mean changes 
in the U and / parameters are consistent with 0. Given 
these results, we use the covariance matrix obtained by 
bootstrapping to characterize the systematic uncertainty 
associated with excluding the p(1700). 
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B. B Background Branching Fractions 

We account for uncertainties on the branching frac- 
tions for the various B background classes by performing 
fits to data after increasing and decreasing the expected 
number of B background events by 10%. The value of 
10% is an estimate of the degree of uncertainty on the 
expected number of events in the B background classes 
with the largest contributions (see Table IT]). We calculate 
a systematic covariance matrix for the U and / parame- 
ters from these fits. Element (j,j) is given by 



^'^ "22' 



(34) 



where 6i is the difference between the value of the ith U 
or / parameter in the fit with increased B background 
contributions, and the fit with decreased contributions. 
The systematic errors (cTgyst = \fGii) associated with the 
B-background branching fractions are given in Table [TV] 
The first column of numbers contains the systematic er- 
rors calculated from the square root of the diagonal ele- 
ments of the covariance matrix while the second column 
of numbers contains the ratio of these uncertainties to 
the statistical uncertainties from the nominal fit to the 
full dataset. 



C. p Lineshapes 

The systematic uncertainties associated with the 
p(770)''', p(770)°, and p(1450) hneshapes are calculated 
by varying their masses and widths according to the un- 
certainties listed in Table Fvl and the correlation matrix in 
Table IVII These correlations and uncertainties were de- 
termined using a Gounaris-Sakurai model fit to data from 
r+ — >■ i7t-/?+(— >■ tt+tt") decays and e+e~ — >■ p° — )■ tt+tt" 
annihilation [15]. We use updated lineshape fits includ- 
ing new data from e+e^ annihilation |16j and r spectral 
functions JL7,. 

The correlations between the three p masses and 
widths as well as their uncertainties were used to cal- 
culate the corresponding covariance matrix. A six- 
dimensional correlated Gaussian was defined with means 
corresponding to the central values of the lineshape pa- 
rameters. A set of 50 vectors of lineshape parameters was 
sampled randomly from the multi-dimensional Gaussian 
distribution and used to perform 50 fits. Using the ran- 
domly sampled masses and widths, and initializing all 
other parameter values to those from the best fit to the 
full on-resonance dataset with nominal parameters, we 
performed a fit with each resonance configuration, pro- 
ducing 50 final sets of \J and / parameters. The covari- 
ances of these 50 sets of V and / parameter values were 
calculated using the set of final fit values from the nom- 
inal fit as expected values. 

As a test of the fitting framework, the mean of the 
number of signal events from the 50 fits was extracted 
and compared to the nominal value. They were found to 



TABLE IV: Square root of the diagonal elements of the covari- 
ance matrix associated with the uncertainties on the U and / 
parameters due to uncertainties on B-background branching 
fractions. The rightmost column contains the ratio of these 
systematic uncertainties to the statistical uncertainties on the 
IJ and 7 parameters as obtained in the nominal fit. 



Parameter 


0"syst 


0"syat/0"stat 


h 


0.0035 


0.09 


/_ 


0.007 


0.11 


rim 


0.012 


0.03 


rRe 


0.017 


0.03 


/+ 


0.0022 


0.04 


rim 
-'+0 


0.025 


0.07 


rRc 
-'+0 


0.031 


0.07 


rim 


0.034 


0.05 


It. 


0.032 


0.04 


Uo 


0.0011 


0.02 


u^ 


0.007 


0.25 


UZo'-^ 


0.14 


0.33 


f/_-6^^ 


0.043 


0.12 


[/+r 


0.07 


0.35 


u+t 


0.031 


0.18 


u- 


0.009 


0.09 


f/+ 


0.0018 


0.03 


rr-,Im 

'-'+0 


0.017 


0.06 


rx-,Re 


0.011 


0.04 




0.016 


0.10 


rj+.Re 

^-,Im 
^-^Rc 

^+,Re 

UT 


0.011 

0.10 

0.019 

0.05 

0.06 

0.012 


0.07 
0.20 
0.04 
0.20 
0.22 
0.13 



TABLE V: p lineshape parameters and uncertainties used in 
the fits and in evaluating the uncertainties on the U and I 
parameters [15i . 



Parameter 


Value [15] (MeV/c^) 


'^p(770)± 


775.5 ± 0.6 


™-p{770)0 


773.1 ±0.5 


rp(770)± 


148.2 ±0.8 


rp{770)° 


148.0 ±0.9 


"lp(1450) 


1409 ± 12 


rp(1450) 


500 ± 37 



be in agreement with a difference of O.So-gtat (where Cstat 
is the statistical error on the number of signal events 
from the nominal fit to the full on-resonance dataset), 
exhibiting negligible bias. As a further comparison, the 
ratios of the systematic errors (taken from the square root 
of the diagonal of the covariance matrix) to the statistical 
errors (taken from the nominal fit) were calculated and 
found to be small as shown in Table IVIll 



17 



TABLE VI: Correlations between p lineshape parameters used 
in evaluating the uncertainties on the U and / parameters [15] . 





n^p{770)± 


"^p(770)0 


rp{770)± 


rp(770)0 


"^p(1450) 


rp(1450) 


'Tlp(770)± 


1.000 


0.109 


0.315 


-0.035 


0.017 


-0.150 


"^p(770)0 


0.109 


1.000 


0.049 


0.290 


0.142 


-0.065 


rp(770)± 


0.315 


0.049 


1.000 


0.361 


0.133 


0.024 


rp{770)0 


-0.035 


0.290 


0.361 


1.000 


0.180 


0.083 


'Tlp(1450) 


0.017 


0.142 


0.133 


0.180 


1.000 


0.779 


rp(1450) 


-0.150 


-0.065 


0.024 


0.083 


0.779 


1.000 



TABLE VII: Square root of the diagonal elements of the co- 
variance matrix associated with the uncertainties on the U 
and I parameters due to uncertainties on the p lineshape. 
The rightmost column contains the ratio of these systematic 
uncertainties to the statistical uncertainties on the U and / 
parameters as obtained in the nominal fit. 



Parameter 


Csyst 


^syst/Cstat 


lo 


0.0005 


0.01 


/_ 


0.0013 


0.02 


rlnx 
-'-0 


0.009 


0.02 


rRe 
-'-0 


0.025 


0.04 


1+ 


0.0007 


0.01 


I'ro 


0.012 


0.03 


rRe 

-'+0 


0.06 


0.13 


rim 
-'-1— 


0.05 


0.07 


7^ 


0.07 


0.09 


Uo 


0.0014 


0.03 


u^ 


0.0020 


0.07 


uzt 


0.029 


0.07 


u-_t 


0.018 


0.05 


f/+r 


0.020 


0.10 


ult 


0.017 


0.10 


u- 


0.0014 


0.01 


u+ 


0.0017 


0.02 


TT-,lin 


0.019 


0.07 


Tj-,nc 


0.011 


0.03 


rr+,Im 


0.007 


0.04 


7-r+,Re 

'-'+0 


0.005 


0.03 


^-,lm 


0.028 


0.06 


^-3c 


0.023 


0.05 


^+,Im 


0.020 


0.08 


^+,Ro 


0.026 


0.10 


u+ 


0.0013 


0.02 



TABLE VIII: Magnitude of changes in U and / parameters 
between fits with and without the uniform background com- 
ponent, the ratio of these magnitudes relative to statistical 
uncertainties from the fit without the uniform component, 
the bootstrapped estimate of the uncertainty on changes in U 
and / parameters between fits with and without the uniform 
component, and the ratio of the mean change in the U and 
I parameters across all bootstrapped fits to their bootstrap- 
estimated uncertainties. 



Parameter 


\AU\ 


lAC/|/<73tat 




(AC/) /a% 


/o 


0.005 


0.13 


0.006 


-0.62 


7_ 


0.006 


0.10 


0.007 


0.55 


7-1 m 


0.041 


0.10 


0.16 


-0.18 


rRe 


0.09 


0.16 


0.25 


0.34 


1+ 


0.0013 


0.02 


0.006 


0.22 


A-0 


0.016 


0.05 


0.08 


-0.11 


rRe 


0.029 


0.06 


0.11 


-0.07 


rim. 


0.006 


0.01 


0.17 


-0.04 


7?1 


0.07 


0.09 


0.21 


-0.15 


Uo 


0.000005 


0.00 


0.008 


-0.07 


K 


0.0036 


0.12 


0.0036 


0.79 


u-_t 


0.14 


0.33 


0.18 


-0.58 


uzi"^ 


0.019 


0.05 


0.13 


-0.03 


t/+6'-" 


0.06 


0.31 


0.06 


0.64 


ult 


0.026 


0.16 


0.06 


-0.33 


u- 


0.006 


0.06 


0.012 


0.39 


ut 


0.005 


0.07 


0.010 


0.30 


TT-,lm 
^+0 


0.0032 


0.01 


0.10 


0.16 


rx-.Re 


0.031 


0.10 


0.11 


0.27 


rr+,Im 
^+0 


0.0026 


0.02 


0.041 


-0.08 


rr+,Rc 


0.006 


0.04 


0.044 


-0.14 


j^-,Im 


0.045 


0.09 


0.16 


-0.38 


^-3c 


0.0032 


0.01 


0.18 


-0.05 


f;+,Im 


0.043 


0.16 


0.06 


-0.34 


jj+,nc 


0.05 


0.21 


0.13 


-0.23 


u+ 


0.007 


0.07 


0.012 


-0.47 



fits to bootstrapped samples are not significant relative 
to the bootstrap-estimated uncertainties (see last column 
of Tabic |VIli| . Also, as in the case of the p(1700) sys- 
tematic, the changes in the U and / parameters with and 
without the extra component are not significant relative 
to the statistical uncertainties on the U and / parameters 



(see middle column in Table VIII I 



D. Uniform Background Contributions 

The systematic uncertainties associated with incoher- 
ent uniform background contributions in the DP are cal- 
culated in the same manner as for the p(1700). The 
nominal fit with the uniform background component al- 
locates only 69.4 events to the uniform component while 
the number of signal events decreases by 28.8 events or 
0.3cr. The changes in the U and / parameters across the 



E. Other Contributions 

In the 2007 BaBar analysis [2], several uncertainties 
were considered and found to provide only small contri- 
butions to the systematic uncertainty: 

• Uncertainties in Arrid, the 5° lifetime, each At 
resolution parameter, tagging fractions, self-cross- 
feed fractions, i3-background tagging fractions, and 
other minor systematics. 
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TABLE IX: Systematic uncertainties for the three sources of 
systematic uncertainty not revisited for the present analysis. 
These values are taken from the 2007 BaBar analysis and used 
without modification [2]. 



Parameter 


B Bkgnd CP Tag 


Side Interference 


"Others" 


lo 


0.004 


0.000 


0.002 


I- 


0.013 


0.003 


0.007 


i-'o 


0.064 


0.018 


0.077 


rRe 


0.083 


0.006 


0.059 


1+ 


0.009 


0.011 


0.005 


I'^ 


0.050 


0.003 


0.065 


rRe 

-'+0 


0.121 


0.019 


0.092 


rim 


0.168 


0.033 


0.133 


jn._ 


0.088 


0.026 


0.078 


Uo 


0.015 


0.000 


0.004 


u^ 


0.005 


0.001 


0.004 


uzi^ 


0.052 


0.007 


0.016 


u-_t 


0.044 


0.022 


0.046 


u+t^ 


0.038 


0.022 


0.011 


rr+.Rc 


0.015 


0.012 


0.007 


u- 


0.041 


0.004 


0.015 


ul 


0.014 


0.003 


0.010 


rr-,Im 

'-'+0 


0.073 


0.028 


0.038 


rr-,R.c 


0.052 


0.004 


0.037 


rr+,Im 
"^+0 


0.042 


0.001 


0.032 


fr+,Rc 


0.059 


0.031 


0.060 


^-,Im 


0.055 


0.031 


0.045 


^-,Ro 


0.238 


0.044 


0.112 


^ + ,Im 


0.028 


0.031 


0.012 


jj + ,Kc 


0.038 


0.028 


0.079 


u+ 


0.036 


0.007 


0.009 



• Uncertainties on CP violation in the B backgrounds 
(calculated by varying the parameters describing 
CP violation in each B background class according 
to their uncertainties). 

• Uncertainties caused by interference between h — > 
cud and b — >■ ucd on the tag side. 

These uncertainties are not expected to provide signif- 
icantly different contributions in the current analysis. 
Therefore, these studies were not repeated, though their 
contributions to the systematic covariance matrix are in- 
cluded. The systematic uncertainties for these contribu- 
tions, as calculated in the 2007 BaBms. analysis, are given 
in Table lIXl 



TABLE X: The square root of the diagonal elements of the 
final systematic covariance matrix for the U and I parame- 
ters, the ratio of these total systematic errors to the statistical 
errors from the final fit, and the ratio of the total error (includ- 
ing statistical and systematic contributions) to the statistical 



Parameter 


O'syst 


CJ'syst/CJ'stat 


\/^s'tat 


+ O'syst/'^stat 


/o 


0.022 


0.59 




1.16 


7_ 


0.028 


0.45 




1.10 


rim 


0.5 


1.08 




1.47 


1% 


0.8 


1.34 




1.67 


1+ 


0.029 


0.46 




1.10 


rim 

-'+0 


0.43 


1.21 




1.57 


rRo 

-'+0 


0.7 


1.57 




1.86 


rim 

'' + - 


0.9 


1.26 




1.61 


7?1 


1.0 


1.30 




1.64 


t^o" 


0.034 


0.65 




1.19 


Uo^ 


0.020 


0.68 




1.21 


uzi"^ 


0.5 


1.19 




1.56 


uzf 


0.37 


1.05 




1.45 


ult 


0.24 


1.16 




1.53 


utt 


0.18 


1.10 




1.49 


u- 


0.06 


0.59 




1.16 


u+ 


0.05 


0.72 




1.23 




0.5 


1.59 




1.88 


rr-,Re 


0.38 


1.19 




1.55 


rr+,Im 


0.17 


1.09 




1.48 


rr+,Ro 


0.19 


1.26 




1.61 


^-,Im 


0.8 


1.68 




1.95 


^-3e 


0.8 


1.70 




1.98 


jj+,lm 


0.26 


1.02 




1.43 


jj+,n. 


0.33 


1.29 




1.63 


u+ 


0.07 


0.77 




1.26 



ratio of the total systematic error to the statistical error 
from the final fit and the ratio of the total error (calcu- 
lated by adding the systematic and statistical errors in 
quadrature) to the statistical error. 

Tables |XI| - |XIII contain the complete systematic corre- 
lation matrix for the U and / parameters. Because the 
three additional sources of systematic uncertainty taken 
from our previous analysis have small contributions, and 
only the diagonal elements of their covariance matrix are 
available to us, we only use the diagonal elements when 
creating the total systematic covariance matrix that is 
used, in turn, to generate the total correlation matrix. 



F. Total Systematic Uncertainties 



VIII. RESULTS 



The square root of the diagonal elements of the total 
systematic covariance matrix for the U and / parameters 
are provided in Table |XJ These values include all system- 
atic uncertainties described above, including the three 
sources of systematic uncertainty that were not recalcu- 
lated in the present analysis. Table Ix] also contains the 



The final values of the U and / parameters extracted 
from the extended maximum likelihood fit to the full 



on-resonance dataset are provided in Table XIV Ta- 
bles |XVf]XVlT| present the statistical correlation matrix 
for the U and / parameters in the fit. From an on- 
resonance dataset containing 53,084 events, the fit ex- 
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TABLE XI: Systematic correlation matrix for the U and / parameters, including contributions from three sources of minor 
systematic uncertainties evaluated in the 2007 BaBar. analysis. Elements above the diagonal are redundant and not included 



in the table. (Continued in Tables XII and XIII 





h 


I- 


rim 
^-0 


rRc 
^-0 


/+ 


rim 

-'+0 


rRc 

^+0 


jto 


/o 


1.00 
















/_ 


-0.17 


1.00 














rim 


0.13 


-0.13 


1.00 












rRc 


-0.07 


-0.00 


-0.03 


1.00 










/+ 


-0.20 


-0.13 


-0.11 


-0.01 


1.00 








rim 
-'+0 


0.04 


-0.09 


0.06 


0.02 


0.00 


1.00 






rRc 

-'+0 


0.21 


-0.05 


0.04 


0.08 


-0.14 


0.16 


1.00 




rim 


0.13 


-0.00 


0.09 


0.09 


-0.08 


0.15 


0.26 


1.00 


1^1 


0.06 


0.08 


0.20 


0.01 


-0.31 


-0.07 


0.16 


0.20 


u,-; 


0.09 


0.04 


-0.13 


0.08 


0.10 


-0.04 


-0.09 


-0.08 


u^ 


-0.11 


0.10 


-0.07 


0.08 


0.09 


-0.08 


-0.05 


-0.12 


f/Io^- 


0.08 


-0.21 


0.10 


-0.31 


0.01 


0.05 


0.01 


-0.11 


f/_-^^° 


-0.03 


-0.01 


-0.41 


0.13 


0.09 


-0.01 


0.03 


-0.16 


c/lo'- 


-0.18 


0.21 


-0.04 


0.18 


0.02 


-0.17 


-0.11 


-0.00 


rx + .Ro 

"-'-0 


0.19 


-0.15 


0.32 


-0.29 


-0.09 


0.00 


-0.03 


0.15 


u- 


-0.12 


0.04 


-0.08 


0.15 


0.05 


-0.02 


0.11 


-0.05 


u+ 


0.01 


-0.08 


-0.23 


0.06 


0.24 


0.27 


0.09 


-0.06 


rr-.lm 


-0.00 


0.07 


-0.11 


0.21 


0.01 


0.04 


0.45 


0.14 


rr-.Rc 


-0.07 


0.11 


-0.12 


0.12 


-0.09 


-0.17 


0.13 


-0.11 


rr+.Im 


0.02 


0.04 


-0.08 


0.11 


-0.02 


0.01 


0.18 


0.09 


fr+,Re 


0.06 


-0.01 


0.08 


-0.02 


-0.07 


0.11 


-0.07 


0.08 


^-.Im 


0.04 


-0.01 


0.19 


0.04 


-0.13 


-0.28 


0.02 


0.19 


^-,Re 


-0.10 


0.07 


0.11 


-0.04 


-0.15 


-0.25 


-0.37 


-0.30 


^+,Im 


0.05 


-0.08 


0.11 


-0.20 


-0.07 


0.08 


-0.22 


0.05 


^+,Hc 


0.24 


-0.11 


0.33 


-0.26 


-0.17 


0.15 


0.07 


0.16 


u+ 


0.09 


-0.04 


0.19 


-0.07 


-0.23 


-0.15 


-0.25 


0.03 





















TABLE XII: (Continued from Table XI I Systematic correlation matrix for the U and / parameters, including contributions 



from three sources of minor systematic uncertainties evaluated in the 2007 BaBar, analysis. Elements above the diagonal are 



redundant and not included in the table. (Continued in Table XIII 





ir- 


Uo 


u^ 


rr-Jm 


t/r„^° 


t/+r 


fr+,Rc 


uz 


ut 


,Hc_ 


1.00 


















u^ 


-0.22 


1.00 
















u^ 


0.07 


0.07 


1.00 














u:t' 


0.05 


-0.14 


-0.17 


1.00 












uzr 


-0.18 


0.11 


-0.04 


-0.09 


1.00 










c/+r 


0.02 


0.08 


0.35 


-0.49 


0.04 


1.00 








f,+.Rc 


0.36 


-0.15 


-0.03 


0.26 


-0.34 


-0.14 


1.00 






u- 


0.01 


-0.07 


0.04 


-0.05 


0.07 


0.05 


-0.22 


1.00 




u+ 


-0.58 


0.24 


0.01 


-0.13 


0.24 


-0.04 


-0.40 


0.05 


1.00 


TT-Mn 

^+0 


0.05 


-0.14 


0.02 


-0.12 


0.09 


0.03 


-0.20 


0.18 


0.12 


r,-,Rc 


-0.07 


-0.04 


0.03 


-0.12 


0.12 


0.10 


-0.27 


0.16 


0.04 




-0.02 


0.11 


-0.18 


-0.12 


0.16 


-0.00 


-0.07 


0.02 


-0.00 


f/+,Rc 


-0.07 


-0.12 


-0.26 


0.06 


-0.12 


-0.21 


-0.01 


-0.02 


-0.02 


j^-,Im 


0.45 


-0.14 


0.06 


0.08 


-0.18 


0.08 


0.35 


-0.13 


-0.57 


^-^Hc 


0.27 


-0.06 


0.16 


0.08 


-0.11 


0.12 


0.17 


-0.06 


-0.48 


U + ,lm 


0.02 


0.02 


-0.26 


0.21 


-0.16 


-0.27 


0.23 


-0.20 


-0.17 


j^ + .Re 


0.33 


-0.15 


-0.13 


0.34 


-0.34 


-0.29 


0.48 


-0.19 


-0.29 


u+ 


0.28 


-0.14 


-0.06 


0.20 


-0.20 


-0.02 


0.30 


-0.23 


-0.42 
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TABLE XIII: (Continued from Table XII I Systematic correlation matrix for the U and I parameters, including contributions 
from three sources of minor systematic uncertainties evaluated in the 2007 BaBar analysis. Elements above the diagonal are 
redundant and not included in the table. 







r;-,Rc 


rr+,Im 
^+0 


rr+,Ro 


y-^Im 


j^-^Ro 


y + Jm 


f^+3c 


u+ 


TT-,i-m 
rr-,Re 


1.00 
0.15 


1.00 
















rr+, Im 


0.06 


-0.01 


1.00 














'-'+0 


-0.04 


-0.03 


-0.10 


1.00 












j^-,Im 


-0.08 


-0.08 


0.07 


-0.04 


1.00 










j^-,Kc 


-0.25 


0.05 


-0.05 


-0.07 


0.18 


1.00 








f; + .lm 


-0.29 


-0.19 


-0.01 


0.05 


0.02 


0.15 


1.00 






^ + .Ho 


-0.16 


-0.21 


-0.13 


0.14 


0.16 


0.09 


0.28 


1.00 




u+ 


-0.19 


-0.16 


-0.05 


0.01 


0.42 


0.28 


0.24 


0.28 


1.00 























tracts 2,940±100 signal events and 46,750±220 contin- 
uum events. The goodness of fit is illustrated in Fig. |6j 
which shows overlaid distributions of fit variables in the 
data used in the final fit and in a parameterized MC sam- 
ple generated using the results of the final fit and equiv- 
alent to 10 times the integrated luminosity of the data 
sample. The signal component of these plots is enhanced 
by a restrictive selection criterion on the NN variable. 



A study of the U and / parameters (see Appendix A 1 ) 
establishes that there is negligible bias in their extrac- 
tion and good robustness in the presence of statistical 
fluctuations. 



A. Quasi-Tw^o-Body Parameters 

The U and / parameters and associated correlations 
can be used to extract the values of the Q2B B^{B ) ^■ 
p^TT^ GP-violation parameters in the time- dependent de- 
cay rate defined in Ref. [7]: 



a±^T 



r^z (At)^(i±^p.) 



-|At|/r 

"4^ 



where Q 



X [1 + Qtag(5 ± A5) sin(AmrfAi) 
-Qtag(C ± AC) cos(Am<iAf)] , (35) 

tag = -t-1 (—1) when the tag-side B meson is a 
B^ {B ). The time- and flavor-integrated charge asym- 
metry ApT^ quantifies direct CP violation, while S and 
C parameterize mixing-induced CP violation related to 
the angle a, and fiavor-dependent direct CP violation, 
respectively. The parameter AC describes the asymme- 

+ T(lf -^ p~Ti+) 
while A5 relates 



try between the rates T{B^ 
and T{B° -^ p--K+) + T{lf 



P+TT 



to the strong-phase difference between the different am- 
plitudes involved in the decay i?° — ^ pir. The U and / 
parameters are related to the Q2B parameters through 
the relations of Eqs. ([l7])-([22]). 

We can also use the U and / parameters and associ- 

p"^n 



TABLE XIV: U and I parameter values from the final fit to 
the complete on-resonance dataset, where the first error is 
statistical and the second is systematic. 



Parameter 


Final Fit Value 


/o 


-0.042 ±0.038 ±0.022 


I- 


-0.00 ±0.06 ±0.03 


/i'S 


-0.61 ±0.43 ±0.46 


rRe 

-'-0 


0.4 ±0.6 ±0.8 


1+ 


0.05 ±0.06 ±0.03 


4" 


-0.04 ±0.36 ±0.43 


rRe 

-'+0 


0.5 ±0.5 ±0.7 


T-Im 


-0.5 ±0.7 ±0.9 


7?1 


-0.6 ±0.8 ±1.0 


u^ 


0.04 ±0.05 ±0.03 


ut 


0.225 ±0.030 ±0.020 


uzt' 


0.53 ±0.44 ±0.52 


f/Io^^ 


0.49 ±0.35 ±0.37 


f/+r 


-0.39 ±0.20 ±0.24 


j^+3c 


-0.05 ±0.17 ±0.18 


u- 


-0.27 ±0.10 ±0.06 


u+ 


1.22 ±0.07 ±0.05 


u+il^ 


0.10 ±0.29 ±0.45 


TT-,-Rc 


0.30 ±0.32 ±0.38 


/■r + .Im 


0.41 ±0.16 ±0.17 


rr+,Rc 


0.01 ±0.15 ±0.19 


U+'"^ 


1.1 ±0.5 ±0.8 


^-^R,c 


-0.5 ±0.5 ±0.8 


f; + ,Im 


-0.07 ±0.26 ±0.26 


^ + ,Re 


-0.19 ±0.25 ±0.33 


u+ 


0.25 ±0.09 ±0.07 



ated correlations to extract the i?° ~> o^tt" CP-violation 



TABLE XV: Statistical correlation matrix for the U and / parameters, 
included in the table. (Continued in Tables |X Vl] and |X VIl[ ) 
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Elements above the diagonal are redundant and not 





h 


I- 


T-Im 
^-0 


rRc 


/+ 


rim 

-'+0 


rRc 

-'+0 


Jim 


/() 


1.00 
















/_ 


-0.04 


1.00 














rim 
-'-0 


0.00 


-0.06 


1.00 












rRc 

-f-o 


-0.05 


0.05 


0.05 


1.00 










1+ 


-0.02 


-0.01 


0.00 


-0.00 


1.00 








rim 
-'+0 


-0.03 


-0.01 


0.01 


-0.01 


0.05 


1.00 






rRc 
-'+0 


-0.12 


-0.02 


0.02 


-0.01 


-0.09 


0.17 


1.00 




rim 


0.01 


-0.06 


0.01 


-0.01 


-0.04 


0.01 


0.01 


1.00 


n- 


0.02 


-0.08 


0.02 


-0.04 


-0.01 


0.03 


0.06 


0.04 


u^ 


-0.04 


0.01 


-0.02 


0.10 


0.00 


-0.05 


-0.04 


-0.00 


K 


-0.08 


0.01 


-0.02 


0.01 


-0.00 


0.03 


0.02 


-0.00 


uzi-^ 


-0.05 


-0.06 


0.01 


0.06 


-0.01 


0.05 


0.09 


0.01 


uz,^^ 


-0.03 


0.02 


-0.17 


-0.01 


-0.00 


0.00 


0.00 


-0.01 


fz+o'-" 


0.00 


0.03 


0.05 


-0.05 


0.01 


-0.05 


-0.09 


-0.01 


j,+ ,He 


0.03 


-0.07 


0.17 


-0.11 


-0.00 


0.03 


0.04 


0.01 


u- 


-0.00 


-0.02 


0.03 


0.02 


0.00 


-0.01 


-0.02 


-0.00 


ut 


-0.03 


0.02 


-0.07 


0.06 


-0.01 


0.01 


-0.01 


0.00 


rr-Am 
^+0 


0.01 


-0.00 


-0.00 


0.01 


0.01 


0.05 


0.26 


0.00 


rr-,Re 


-0.08 


0.01 


-0.01 


0.01 


-0.06 


-0.06 


-0.02 


-0.00 


rr+,Im 


-0.01 


-0.01 


-0.00 


0.01 


-0.05 


0.06 


-0.01 


-0.00 


r/+,Re 


0.05 


-0.00 


-0.00 


0.01 


-0.03 


0.09 


-0.09 


-0.00 


^-,Im 


0.01 


-0.04 


0.00 


0.01 


-0.04 


0.01 


0.03 


0.04 


^-,Re 


0.01 


-0.01 


0.00 


-0.01 


-0.00 


0.01 


0.01 


-0.11 


^+,Im 


0.02 


-0.02 


0.02 


-0.03 


-0.01 


0.04 


0.06 


-0.04 


^+,Hc 


0.02 


-0.04 


0.02 


-0.04 


-0.03 


0.04 


0.08 


0.08 


u+ 


0.01 


-0.01 


0.00 


-0.01 


-0.05 


-0.03 


-0.04 


-0.00 





















TABLE XVL (Continued from Table |XV[) Statistical correlation matrix for the U and / parameters. Elements above the 



diagonal are redundant and not included in the table. (Continued in Table XVII 





Ifl 


Uo 


u^ 


UZi'-^ 


uzr 


t/+r 


rr+,Rc 


uz 


u+ 


jRc 


1.00 


















u^ 


-0.03 


1.00 
















u^ 


-0.03 


0.07 


1.00 














uzil"^ 


0.08 


0.02 


0.00 


1.00 












uzr 


-0.02 


-0.05 


0.02 


0.09 


1.00 










[/+o'" 


-0.09 


0.02 


0.05 


-0.26 


-0.04 


1.00 








u+t 


0.06 


-0.04 


-0.07 


0.06 


0.01 


-0.03 


1.00 






u- 


0.01 


0.00 


-0.03 


-0.08 


-0.07 


0.06 


-0.02 


1.00 




u+ 


-0.05 


0.05 


0.26 


0.02 


0.05 


-0.05 


-0.05 


-0.12 


1.00 


Tr-,l-ni 


0.01 


0.13 


0.05 


0.02 


-0.01 


-0.01 


-0.00 


-0.01 


0.01 


rr-,Re 
"^+0 


-0.02 


-0.10 


0.04 


-0.03 


0.02 


0.03 


-0.02 


0.01 


0.03 


rr+,Im 


-0.01 


0.04 


0.15 


0.03 


0.02 


-0.01 


-0.00 


-0.01 


0.07 


r, + ,Rc 


-0.01 


0.02 


-0.09 


-0.00 


0.01 


-0.00 


0.01 


-0.00 


0.02 


(^-,Im 


0.05 


0.00 


0.02 


0.06 


0.03 


-0.05 


0.01 


0.01 


0.06 


^-^Re 


-0.04 


-0.01 


-0.01 


0.02 


0.00 


-0.02 


0.01 


-0.08 


-0.04 


[; + .Im 


0.01 


-0.03 


-0.01 


0.12 


0.01 


-0.12 


0.06 


0.01 


0.02 


jj+,H. 


0.14 


-0.04 


-0.03 


0.11 


-0.01 


-0.12 


0.07 


-0.07 


-0.06 


u+ 


0.04 


-0.01 


0.01 


0.04 


0.01 


-0.04 


0.01 


-0.05 


0.04 
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m(ii* It) (GeV/c') 




3 4 5 
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3 4 5 

m(;t-ji°)(GeV/c') 



FIG. 6: (color online) Overlay plots of fit variable distributions in on-resonance data (points with error bars) and in a pa- 
rameterized MC sample generated from the final fit results (red line) with 10 times the number of events in data. The MC 
histograms are scaled to have the same integral as the data histograms onto which they are overlaid. A restrictive selection 
criterion is applied to the NN output (NNncw > 0.8) to enhance the signal component. 



parameters and decay fraction: 



C/n 



r — 



^00 — 



r7+' 

2A, 

f/+' 



■^°" U+ + U1 + u+ ■ 



(36) 

(37) 
(38) 



These eight Q2B parameters related to p^n^ and p^vr^ 
decays are extracted using a x^ minimization technique 
that accounts for the statistical and systematic correla- 
tions between the U and / parameters. In each step of 
the minimization process, the current values of the free 
parameters are used to calculate the corresponding values 
of the eight U and / parameters on which they depend 



and a vector Vdig is constructed from the differences be- 
tween these U and / parameter values and the values 
obtained in our final fit to the full on-resonance dataset. 
The x^ is then calculated as 



XQ2B 



i^dTff(c^''"*")"'^d 



iff, 



(39) 



where C"^^**^ is the 8x8 covariance matrix for the rel- 
evant U and / parameters from the fit to data. Ta- 



ble XVIII presents the Q2B parameters extracted from 
the full fit to on-resonance data along with their statisti- 
cal and systematic errors. In a study of correlations be- 
tween a At and DP position we find a small contribution 
to the systematic uncertainty on /oo, which is included 
in Table |XVIII[ The correlation matrix for the Q2B pa- 
rameters is provided in Table |XIX[ A study of the Q2B 



parameters (see Appendix A2) establishes that there is 
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TABLE XVII: (Continued from Table XVI I Statistical correlation matrix for the U and I parameters, 
diagonal are redundant and not included in the table. 



Elements above the 







rr-,Rc 




rr+,Ro 


^-Jm 


f^-Jlc 


[; + Jm 


^+Mc 


u^ 


TT-,lm 


1.00 


















rr-,Re 


-0.14 


1.00 
















rr+,Im 


0.22 
0.03 


0.07 
-0.01 


1.00 
-0.09 


1.00 












j^-,lm 


-0.00 


0.00 


0.03 


0.00 


1.00 










j^-,He 


0.00 


-0.00 


-0.00 


-0.00 


0.31 


1.00 








^+.lm 


0.01 


-0.02 


0.02 


-0.00 


0.18 


-0.06 


1.00 






^ + ,Hc 


0.01 


-0.03 


0.00 


-0.00 


-0.17 


0.06 


0.13 


1.00 




u+ 


0.04 


-0.01 


-0.00 


0.02 


-0.14 


0.01 


0.06 


0.01 


1.00 























TABLE XVIII: Quasi-two-body parameter values and uncer- 
tainties corresponding to the fit to the complete on-resonance 
dataset. 



Par am 


Value CTstat O-syst 


Apn 


-0.100 0.029 0.021 


C 


0.016 0.059 0.036 


AC 


0.234 0.061 0.048 


S 


0.053 0.081 0.034 


AS 


0.054 0.082 0.039 


Coo 


0.19 0.23 0.15 


Soo 


-0.37 0.34 0.20 


foo 


0.092 0.011 0.009 



negligible bias in their extraction and good robustness in 
the presence of statistical fluctuations. 

The parameters ApT^ and C can be transformed into 
the direct CP-violation parameters Ap~ and A~t^ where 



a: 



A-+ 



T(B" 



p n^ 



r(B" 



P^TT 



r(i?° 


-^ 


p 


n+) 


+ r(BO 


->p+ 


n-) 


t{b' 


-^ 


P^ 


^-) 


-r(B" 


-^ p- 


n+) 



T{B 
using the relations 



p+Tr-)+T(BO ^ p-TT+)' 



A+- = - 



^p7T 



A-+ - 



^pTT 



1 + AC + Ap^C ' 
1-AC-Ap^C 



(40) 
(41) 

(42) 
(43) 



We extract the central values and uncertainties for 
these parameters using a x^ minimization procedure in 
the two-dimensional plane corresponding to A'^~ versus 
A~^. At each point in the plane, the values of A'j^^ and 
A~^ are fixed and used in combination with AC (which 
is free to vary) to determine the corresponding values of 
Ap^ and C. These values are then used in combination 
with AC and the five other Q2B parameters to calculate a 



X^ value as described above. From this two-dimensional 
scan, we find 



A 



pTT 



0.09+[!;J]^ ± 0.04, 



^;+ = -o.i2±o.08trot 



(44) 
(45) 



with a correlation of 0.55 evaluated from the la con- 
tour for statistical and systematic errors combined. A 
two-dimensional likelihood scan with combined statisti- 
cal and systematic uncertainties and 68.3%, 95.5%, and 
99.7% confidence-level contours is provided in Fig.JTJ The 
origin, corresponding to no direct CP violation, lies on 
the 96.0% confidence-level contour (Ax^ — 6.42), corre- 
sponding to a p = 4.0% probability, under the assump- 
tion of no direct CP violation, of obtaining a result that 
deviates from the origin at least as much as ours. 



B. Alpha Scan 

In order to extract likely values of a from the U and / 
parameters obtained in our final fit, we perform a scan of 
a from 0° to 180°. At each scan point, a x^-minimization 
fit is performed using the goodness of fit measure: 



x'(«) = [v 



data 



_ Tyscaii -* 

rr+ scan 



(C 

, 2 



date 



)-i [V 



data 



y« 



ml 

(46) 



where C"^^*'' is the 26 x 26 covariance matrix for Us and 
Is from our fit to data and y^***^ and V^^^'^ are vectors of 
the 26 U and / parameters from our fit to data and the 
current iteration of the minimization, respectively. The 
last term is a Gaussian constraint that restricts U^ in 
the scan to lie within e of 1 (the value to which it is fixed 
in the fits to data); e is set to 0.0001. Because the overall 
scaling of the U and / parameters is not physically mean- 
ingful, U^ is allowed to be exactly 1 in the fit and this 
constraint term does not significantly contribute to the 
total x^ ■ III the x^-minimizations, the actual free param- 
eters are the tree (T*''^) and penguin (P^) amplitudes. 
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TABLE XIX: Combined statistical and systematic correlation matrix for the quasi-two-body parameters corresponding to the 
fit to the complete on-resonance dataset. Values above the diagonal are redundant and omitted for clarity. 





v^pTT 


C 


AC 


S 


AS 


Coo 


5oo 


/oo 


Apn 


1.000 
















c 


0.035 


1.000 














AC 


0.154 


0.213 


1.000 












S 


-0.040 


-0.065 


-0.070 


1.000 










AS 


-0.041 


-0.038 


-0.060 


0.199 


1.000 








Coo 


-0.088 


-0.041 


-0.034 


0.026 


0.011 


1.000 






Soo 


-0.005 


0.007 


0.044 


-0.081 


-0.007 


0.002 


1.000 




/oo 


0.074 


0.009 


0.016 


0.029 


-0.016 


-0.062 


0.062 


1.000 





















+ Cl. 




where 



-0.75 -0.50 -0.25 



0.25 0.50 0.75 



4" + 



FIG. 7: (color online) Combined statistical and systematic 
two-dimensional likelihood scan of .4p.^ versus Ap^ with 
68.3%, 95.5%, and 99.7% confidence-level contours [Ax^ = 
{2.30,6.18,11.83}). The yellow dot inside the contours indi- 
cates the central value. 



which are related to the p-resonance amplitudes through 
the formulae 



A+ 


= T+e-"' + P+, 


(47) 


A- 


= T-e-"'+p-, 


(48) 


A" 


= T°e-*" + P°, 


(49) 


a' 


= T-e+'" + p-, 


(50) 


A^ 


^ T+e+'"+P+, 


(51) 


a" 


= T°e+''" + P", 


(52) 



A+ 


EE A{B^ ^ p+vr-), 


(53) 


A- 


= A(B" ^ p-7r+). 


(54) 


A° 


= A(i3" ^ pVO), 


(55) 


3+ 


P 


(56) 


A^ 


^ ^-A(B"^p-n+), 
P 


(57) 


a' 


P 


(58) 



Note that due to SU{2) flavor symmetry, the third pen- 
guin amplitude can be calculated from the other two us- 
ing the relation P" = -i (P+ + P") (see Refs. [H] and 
[20]). At each step in the minimization process, the cur- 
rent values of the tree and penguin amplitudes as well as 
the current fixed value of a are used to calculate the p- 
resonance amplitudes. These p amplitudes are then used 
to calculate the U and / parameters that comprise the 
vector y'^^.n^ using Eqs. (12 1 through (16), which relate 



the U and / parameters to the p amplitudes. In the fits, 
we take advantage of a global phase that is not physically 
meaningful to fix the phase of T"*" to 0. 

As the scan proceeds, a minimum x^ value is extracted 
from the fit at each value of a. We convert these x^ values 
to "S" values by calculating the x^ probability of each 
value according to 



f{x; l)dx, 



(59) 



where a is the difference between the x^ ^t the current 
scan point and the minimum x^ for all the scan points, 
and /(x; 1) is a x^ distribution with one degree of free- 
dom. The variable "S" corresponds to what is commonly 
referred to as "1— Confidence Level" (1— C.L.) and is sim- 
ply the p- value of a x^ test at each scan point. 



1. Incorporating Information from Charged B Decays 

Following the methods employed in Belle's 2007 P° — )■ 
pTT analysis |3| and described in Ref. [B], we perform a 
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further a scan that makes use of measurements from the 
charged decays B^ — )■ p*'°7r°'^. Amphtudes for these 
modes can be related to amphtudes in the neutral B 
modes through isospin relations. These relations result 
in four "constraint" equations while introducing only two 
new free parameters in the fit (the real and imaginary 
parts of a tree amplitude, r+°). The charged B measure- 
ments of interest are the branching fractions and asym- 
metries: 



B{p+n^) = c{\A+r + \A-r)rB^, 


(60) 


B{p\+) = c{\A^+\' + \A^-\')rs., 


(61) 


\A-0\2 \A+0\2 
^(P ^ ) - |^-0|2 + |^+0|2' 


(62) 


140-12 140+12 

A(o\+] - ' ' "' ' 


(63) 


where c is a constant and 




A±o = ^A{B^ ^ p±^"), 


(64) 


P 


(65) 



In the fit, we fix c = 1 and no longer require Ut — 1. This 
is equivalent to letting c be a free parameter in the fit and 
setting U^ = 1. Due to this convention, it is necessary to 
divide all the current values of the U and / parameters 
during the minimization process by the current value of 
U^ before using them to calculate the current x^ value. 
According to SU{2) isospin symmetry, several relations 
hold between the amplitudes in Eqs. (53H(58) and (64|- 



(651: 



A+ + A- 



2A" 



2A 



V2{A+''+A"+) 
V2e^^"'{A-" + A^-), 



A+" - A"+ - V2{A+ - A-) 



A- 



A"- 



-V2{A- - A-^ 



(66) 

(67) 
(68) 



(69) 



Based on these relations, we can parameterize the 
charged B amplitudes according to 



V2A+" 
V2A°+ 



e-'"T+° + P+ - P- , 
e"'"(r++T" +2T" -T^ 
-P+ + P-. 



%/2yl"- ^ e+*"(T++r-+2T°-r+°) 
-P+ + P-. 



(70) 

(71) 
(72) 

(73) 



that implicitly incorporates the isospin relations between 
the different charged and neutral B modes. Because the 
global phase is physically irrelevant, we fix it by requiring 
that ImT+ = 0. 

The isospin "constrained" and "unconstrained" scans 



are performed identically, except that Eqs. (70)-(73l are 
included as Gaussian constraints in the x^ calculation. 



The system of equations (60)-(63) is used to express 
the magnitudes of A^^ and A^^ in terms of the branch- 
ing fractions and asymmetries for the charged B modes 
B^ -^ p^n° and B^ — >■ p^n^. Using world average 
measurements from |llj . we calculate the value of each 
of these magnitudes as well as their uncertainties. At 
each step in the minimization process, and for each of 
Eqs. (I70|-(l73l), a term is added to the x^- 



\Ai, 



|A„ 



'^\A„ 



.1 



(74) 



where |v4itcr.| is the magnitude of the relevant A^' or 
A"* parameter for the current iteration of the minimiza- 
tion process, jAmeas. | is the magnitude of the amplitude 
based on branching fractions and asymmetry measure- 
ments, and CTi^^^^^ I is the uncertainty in the value of 
the magnitude due to measurement uncertainties for the 
branching fractions and asymmetries. For those branch- 
ing fractions and asymmetries that have asymmetric un- 
certainties, we choose whether to use the upper or lower 
error for the calculation in a given iteration by ascertain- 
ing whether the value of the branching fraction or asym- 
metry corresponding to the tree and penguin parameters 
in the current iteration of the minimization process is less 
than or greater than the experimental value, respectively. 



2. Results of a Scan 

Plots of the x^ values from our final a scans with 
isospin constraints (solid red) and without isospin con- 
straints (dashed black) are shown in Fig. [SJ The cor- 
responding S distributions are presented in Fig. [9] As 
indicated by our robustness studies (see Appendix), the 
S scan is not robust and cannot be interpreted in terms 
of Gaussian statistics. 



IX. CONCLUSION 

We have performed a time-dependent Dalitz plot anal- 
ysis of the mode i?" — ?► (pTi")" in which we extract 26 U 
and / parameter values describing the physics involved, 
as well as their full correlation matrix in an extended un- 
binned maximum likelihood fit. From these fit results, we 
extract the Q2B parameters v4p^, C, AC, S, AS, Coo, ^oo, 
and /oo, with values given in Table XVIII These Q2B 



values are consistent with the results of the 2007 BaBAR 
analysis [2] as well as the results obtained by the Belle 
the neutral i? amplitudes, this yields a parameterization collaboration [3]. We also perform a two-dimensional 



When combined with Eqs. ( 47 1-( 52 ) , which parameterize 
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" o ' ' ' ' 50 ' ' ' ' 100 ' ' ' ' 150 ' ' « (deg) 

FIG. 8: (color online) Isospin-constrained (solid red) and un- 
constrained (dashed black) scans of minimum x^ values as a 
function of a. The scans are based on the fit to the full on- 
resonance dataset and include contributions from both sta- 
tistical and systematic uncertainties. Note that the origin is 
suppressed. 




- a (deg) 



FIG. 9: (color online) Isospin-constrained (solid red) and un- 
constrained (dashed black) scans of E values as a function of 
a. The scans are based on the fit to the full on-resonance 
dataset and include contributions from both statistical and 
systematic uncertainties. The upper and lower horizontal 
dashed lines correspond to E = 0.05 and 0.32, respectively. 



likelihood scan of the direct CP- violation asymmetry pa- 
rameters for B° —J' p^TT^ decays, finding the change in 
X^ between the minimum and the origin (corresponding 
to no direct CP-violation) to be Ax^ = 6.42. Finally, 
we perform one-dimensional likelihood-scans of the uni- 
tarity angle a (see Figure Isl) both with (solid red) and 
without (dashed black) isospin constraints. However, as 
indicated by our robustness studies (Appendix | A 3[ ), the 
extraction of a with our current sample size is not ro- 
bust. Maximum likelihood estimators are known to be 
Gaussian in general only in the limit of large data sets. 
This analysis would benefit greatly from increased sam- 
ple sizes available at higher-luminosity experiments. 
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Appendix A: Robustness Studies 

We assess the robustness with which the fit frame- 
work extracts statistically accurate values and uncertain- 
ties for the U and / parameters, the Q2B parameters, 
and a, by employing MG-simulated samples generated 
with a parameterized detector simulation and with sig- 
nal and background contributions corresponding to those 
expected in the full on-resonance dataset. The simulated 
samples are generated using physical \J and / parame- 
ters based on specific tree and penguin amplitudes and 
a = 89° (approximately the world average). Each simu- 
lated dataset is generated with the same parameter val- 
ues, but a different random-number seed. By examining 
the results of the fits to each of these simulated datasets, 
we assess the the fit robustness. 

For these studies, we use 25 samples generated with 
different seeds, so that the uncertainty on the bias is 
negligible compared to the other uncertainties. 



1. U and I Parameter Robustness Studies 

For each of the 26 U and / parameters extracted in the 
fits to the MG samples, we calculate the RMS (across the 
25 MC samples) of the differences (measured in units of 
the statistical uncertainty on the extracted values of the 
parameters) between the generated and extracted values 
of the parameters. We also calculate the average across 
the 25 MG samples of these differences for each of the V 
and / parameters. The RMS difference, averaged across 
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TABLE XX: Results of robustness study of quasi-two-body 
parameters. The first column of numbers is the ratio of the 
square root of the variance divided by the mean error on the 
fit to the parameter across 25 MC samples. The second col- 
umn is the average difference (measured in units of statistical 
uncertainty, o, on the extracted values of the parameters) be- 
tween the extracted value of the parameter and its generated 
value, across all MC samples. The third column is the same as 
the second, except that the absolute value of each difference 
is taken before averaging. 



TABLE XXL Results of parameterized MC study of robust- 
ness of a scans. The generated value a is agen = 89° and the 
columns are described in the text. 



Param 


VVariance/ (a) 


Avg #a 


Avg Abs #a 








Diff From 


Diff From 








Gen Val 


Cen Val 


•^pTT 




0.94 


-0.13 


0.76 


c 




1.15 


0.06 


0.90 


AC 




0.94 


0.04 


0.75 


S 




1.11 


0.03 


0.92 


AS 




1.02 


-0.20 


0.82 


Coo 




1.15 


-0.10 


0.89 


5oo 




1.13 


0.23 


0.92 


/oo 




1.08 


0.28 


0.93 



the 26 U and / parameters, is 1.17, while the difference 
in units of statistical uncertainty, averaged across the U 
and / parameters, is 0.04. This demonstrates that the U 
and / parameter values are extracted robustly and with 
negligible bias. 



a Scan 


Upper Lower 
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ij^a Between 


Peak 


Error 


Error 


Error 


x' 
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X^(«gen) 


Gen and Fit 
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43 
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44 


+5 
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-5 


5 
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2.3 


-8.2 


49 


+5 


-5 


5 


24.2 


2.9 


-8.7 


52 


+5 


-5 


5 


16.1 


2.4 


-8.2 


53 


+5 


-21 


13 


23.2 


1.5 


-7.2 


60 


+14 


-8 


11 


16.5 


3.3 


-2.1 


74 


+7 


-6 


7 


21.2 


0.7 


-2.1 


74 


+5 


-13 


9 


15.9 


4.1 


-2.9 


75 


+10 


-8 


9 


21.3 


1.8 


-1.5 


76 


+8 


-18 


13 


21.3 


2.5 


-1.7 


80 


+5 


-6 


6 


30.0 


2.3 


-1.7 


83 


+6 


-7 


7 


24.0 


1.3 


-0.9 


84 


+6 


-6 


6 


26.4 


1.3 


-0.9 


84 


+7 


-6 
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30.1 


1.0 


-0.7 


87 


+7 


-7 


7 


22.9 


0.4 


-0.3 


88 


+7 
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10.9 
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-0.1 


89 
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23.4 


0.00 
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94 
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10.0 


1.3 


0.9 


112 


+5 


-6 


5 


19.0 


3.0 


3.9 


115 


+5 


-6 


5 


23.3 


2.2 


4.6 


124 


+28 


-15 


22 


25.6 


1.4 


2.3 



2. Quasi- Two-Body Robustness Studies 



a Robustness Studies 



The results of the robustness study of the Q2B param- 
eters are provided in Table |XX| In this table, the first 
column of numbers is the ratio of the square root of the 
variance of the parameter across 25 MC samples divided 
by the mean error on the fit to the parameter across all 
MC samples. As one would expect, the square root of the 
variance of each parameter value is approximately equal 
to the mean statistical uncertainty on the variable as ex- 
tracted from the fits to MC samples. The second column 
is the average difference (measured in units of statistical 
uncertainty, a, on the extracted values of the parame- 
ters) between the extracted value of the parameter and 
its generated value, across all MC samples. These values 
are all within 0.28cr of 0, indicating negligible bias. The 
third column is the same as the second, except that the 
absolute value of each difference is taken before averag- 
ing. These values are no greater than 0.93ct, indicating 
that the extracted values of the parameters are reliably 
close to the generated values. Taken together, these re- 
sults indicate that the extraction of the Q2B parameters 
is robust to possible non-Gaussian fluctuations, and un- 
biased. 



The results of the robustness scans of a are provided 
in Table |XXI[ sorted by the absolute difference between 
the extracted a-scan peak position and the generated 
value of 89°. The flrst column lists the position of the 
most favored value of a in the scan. The second and 
third columns list the upper and lower errors, respec- 
tively, which are calculated as the number of degrees to 
either side of the a scan peak position at which the E 
value drops to 0.32. The fourth column lists the mean 
of the upper and lower errors while the fifth column lists 
the value of the minimum x^ obtained at the a scan 
peak position. The second to last column gives a mea- 
sure of the consistency between the likelihoods for the 
peak a position and the generated a position based on 
the change in x^, and the last column is the distance in 
a between the generated and extracted a peak positions 
(where the upper or lower error is used as appropriate). 
For 17 of the 25 scans, the extracted value of a lies within 
3(7 of the generated value. Examining the individual a 
scans reveals three distinct solutions for a that tend to 
be favored (including the generated value of 89°) and 
each scan tends to include at least one secondary peak in 
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addition to the primary peak. Figure 10 illustrates the 
three solutions for alpha by providing the sum of 25 nor- 
malized Gaussians with means and widths determined by 
the peak positions and symmetric errors extracted from 
the 25 a scans. Also plotted are the individual Gaussians 
that contribute to the total PDF. Because the errors are 
not truly Gaussian, Figure [TO] provides an incomplete pic- 
ture of the scan results. A better illustration is provided 
by Figure [TT] which displays the total E distribution ob- 
tained by summing all 25 S a scans after normalizing 
each to the same area. The total distribution is scaled so 
that it peaks at 1. Also plotted are the individual scaled 
S scans that contribute to the total distribution. The fi- 
nal PDF closely resembles that obtained by naively sum- 
ming Gaussian distributions, though it exhibits more fine 
features. Again, the distribution indicates three distinct 
solutions for a, with the generated value of 89° being 
favored. At the la level (I]=0.32), the total scan distri- 
bution allows both the central and left peak. The pres- 
ence of these secondary solutions indicates that with the 
current signal sample size and background levels, there is 
still a significant possibility that the favored value of a in 



a particular scan will correspond to a secondary solution. 

Arbitrary Units 
0.6 




- a (deg) 



FIG. 10: (color online) Plot of the sum of 25 normalized 
Gaussians with means and widths determined by the peak 
positions and mean errors extracted from the 25 a scans of 
parameterized MC generated with a = 89°. Also plotted are 
the individual Gaussians that contribute to the total PDF. 
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